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ABSTRACT 


Jet pumps are devices capable of pumping 
fluids to a higher pressure by inducing the motion of a 
secondary fluid employing a high speed primary fluid. 
The main components of a jet pump are a primary 
nozzle, secondary fluid injectors, a mixing chamber, a 
throat, and a diffuser. The work described in this 
paper models the flow of a two-phase primary fluid 
inducing a secondary liquid (saturated or subcooled) 
injected into the jet pump mixing chamber. The model 
is capable of accounting for phase transformations due 
to compression, expansion, and mixing. The model is 
also capable of incorporating the effects of the 
temperature-and pressure-dependency in the analysis. 
The approach adopted utilizes an isentropic constant- 
pressure mixing in the mixing chamber and at times 
employs iterative techniques to determine the flow 
conditions in the different parts of the jet pump. 

NOMENCLATURE 

Latin Symbols : 

A cross-sectional area, m 2 

a speed of sound, m/s 

h specific enthalpy, kJ/kg 

M Mach number, dimensionless 
m mass flow rate, kg/s 

P pressure, kPa 

s specific entropy, kJ/kg. K 

T temperature, °C 


V velocity, m/s 

x quality 

Greek Symbols : 

$ entrainment ratio, m/rr^, dimensionless 
p density, kg/m 3 

Subscripts : 

d diffuser 

de diffuser exit 

e exit 

f saturated liquid condition 

g saturated vapor condition 

i inlet 

m mixture 

me mixing chamber exit 

mi mixing chamber inlet 

n primary nozzle 

ne primary nozzle exit 

ni primary nozzle inlet 

o stagnation condition 

p primary fluid 

s secondary fluid 

se secondary fluid exit 

si secondary fluid inlet 

ss conditions downstream of shock wave 
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INTRODUCTION 

Jet pumps are devices capable of pumping 
fluids to a higher pressure employing a 

nozzle/diffuser/mixing chamber combination. A 
primary fluid is usually allowed to pass through a 
converging-diverging nozzle where it can accelerate to 
supersonic speeds at the nozzle exit. The relatively 
high kinetic energy that the primary fluid possesses at 
the nozzle exit is accompanied by a low pressure region 
in order to satisfy momentum conservation 

requirements. The low pressure region downstream of 
the nozzle exit permits a secondary fluid to be entrained 
into and mixed with the primary fluid in a mixing 
chamber located downstream of the nozzle. Several 
combinations may exist in terms of the nature of the 
primary and secondary fluids in so far as whether they 
are single- or two-phase fluids. Depending on this, the 
jet pump may be classified as gas/gas, gas/liquid, 
liquid/liquid, two-phase/liquid, or other similar 
combinations. 

The mixing chamber serves to create a 
homogeneous single-phase or two-phase mixture which 
is ready to enter a diffuser where the high kinetic 
energy of the mixture is convened into pressure energy. 
If the fluid mixture entering the diffuser is in the 
supersonic flow regime, a normal shock wave may 
develop inside the diffuser. If the fluid mixture is one 
that can easily undergo phase transformations, a 
condensation shock may develop. Because of the 
overall rise in pressure in the diffuser as well as the 
additional rise in pressure caused by the shock wave, 
condensation of the fluid mixture becomes more likely. 
If the two-phase flow entering the diffuser is 
predominantly gaseous with liquid droplets suspended 
in it, chances are that it would transform into a 
predominantly liquid flow containing gaseous bubbles 
(bubbly flow). 

Because of the high speed of the primary fluid, 
there exists a region known as the potential core region 
where entrainment and mixing between the primary and 
secondary fluids does not take place. Thus, as the 
primary fluid enters the mixing section and expands, 
and because of the lack of mixing in the potential core 
region, the peripheral annular area available for the 
secondary fluid will shrink. Because of continuity 
considerations, the secondary fluid will accelerate in the 
annular region and it would be conceivable that it 
reaches sonic speed there. This secondary fluid 


acceleration will be accompanied by a reduction in 
pressure in the peripheral area which would in turn 
make it more likely for the otherwise subcooled or 
saturated secondary fluid to flash into vapor, thus 
forming either a two-phase mixture or pure vapor in the 
region occupied by the secondary fluid. Now that the 
secondary fluid may have become predominantly or 
purely vapor, any further motion downstream in the 
region occupied by the secondary flow is likely to be in 
the supersonic regime and might be accompanied by 
shock waves in the peripheral region. As the potential 
core region ends, the now two-phase or pure vapor 
secondary fluid will start a mode of mixing with the 
already two-phase primary fluid in the remainder of the 
mixing section downstream of the potential core region 
and will form a homogeneous two-phase mixture which 
is bound towards the diffuser. As stated earlier, this 
two-phase mixture will undergo a velocity reduction 
accompanied by a pressure rise and will be subject to 
undergoing a shock wave somewhere in the diffuser. 

Literature dealing with jet pumps is abundant 
and covers a very wide array of application areas. 
Example application areas includes vacuum pumps 
which are used in the food industry, power station 
work, and the chemical industry; ejector systems which 
have applications in the aircraft industry as cabin 
ventilators and for purposes of jet thrust augmentation; 
jet pumps which are used in the oil industry for oil well 
pumping; and steam-jet ejector refrigeration, to just 
name a few. Examples of work relevant to this 
investigation includes those of Holmes et al. (1987), 
Cunningham and Dopkin (1974), Eiger et al. (1991), 
Marini et al. (1992), Jiao et al. (1990), Neve (1991), 
Bredikhin et al. (1990), and Fairuzov and Bredikhin 
(1995). 

One of the primary motivations behind this 
work is to attempt to understand the performance of a 
two-phase jet pump which constitutes part of a flow 
boiling test facility at NASA-Marshall. The flow 
boiling apparatus is intended to provide data necessary 
to design highly efficient two-phase thermal control 
systems for aerospace applications. The facility will 
also be capable of testing alternative refrigerants and 
evaluate their performance using various heat 
exchangers with enhanced surfaces. The test facility is 
also intended for use in evaluating single-phase 
performance of systems currently using CFC 
refrigerants (Holladay and Hunt, 1996). 



3 


While past researchers have been able to model 
the two-phase flow jet pump using the one-dimensional 
assumption with no shock waves and no phase change* 
there is no research known to the authors apart form 
that of Anand (1992) that was able to account for 
condensation shocks. The work performed by Faimzov 
and Bredikhin (1995), while providing useful insights 
into the performance of the jet pump under study in this 
paper, did not fully account for all the issues involved 
in the constant-pressure mixing process performed in 
the mixing chamber. Thus, one of the objectives of the 
work described here is to model the jet pump in 
question accounting for phase transformations due to 
expansion, compression, and mixing, and to provide a 
more accurate description of the constant-pressure 
mixing process. The analysis employs an isentropic 
homogeneous expansion/compression model, and while 
the results are presented for R-134a, the analysis is 
general enough to be applied to other refrigerants such 
as R-22 and R-l 13. 

ANALYSIS 

The analysis described in this paper deals with 
a jet pump whose primary fluid is a two-phase mixture 
and whose secondary fluid is either a subcooled or 
saturated liquid having the same chemical composition 
as the primary fluid. The primary fluid enters a 
converging-diverging nozzle (primary nozzle) at inlet 
conditions P^, T m , and x^ for the pressure, 
temperature, and quality, respectively, and exits at a 
high velocity and low pressure (see Figure 1). The low 
pressure at the primary nozzle exit induces the flow of 
the secondary liquid (which may also enter by the 
forced action of subcooled liquid spray injectors). The 
secondary liquid is introduced from an upstream 
pressure of P si (and a low velocity) and enters the 
mixing section at a pressure of P^. The constant- 
pressure mixing process as modeled in this paper 
requires the pressure P^ to equal the mixing chamber 
exit pressure P re . This type of mixing has been 
selected since it produces compression ratios in the jet 
pump which are a bit higher than those produced under 
constant-area mixing conditions. 

In order to be able to perform the analysis as 
stated above, the inlet conditions of both the primary 
and secondary fluids as well as the entrainment ratio 
and pump geometry (except the mixing chamber exit 
area) are assumed known. The mixing chamber exit 
area, A w , is assumed unknown and will be determined 


as part of the solution in order to satisfy momentum 
conservation in the mixing chamber, while imposing the 
constant-pressure mixing condition. This situation 
makes the model described here more design (as 
opposed to analysis) oriented. In what follows, a 
summary of the equations employed in the model will 
be presented. 

Primary Nozzle : 

The primary nozzle is segmented into two parts 
for purposes of this analysis; the converging section and 
the diverging section. In both sections, the isentropic 
homogenous expansion model is employed. Since the 
inlet pressure P*, temperature T^, and quality x^, are 
all known, the specific volume, enthalpy, and entropy 
can be determined: 

V m = V ntf + ( V ni* “ 'W U) 


hni = + (*Vg - Kj> (2) 


S ai = S ntf + \i ( S ni^ ~ S ni j) (3) 

Applying the continuity and energy equations 
in the converging section along with the definitions of 
the Mach number and the speed of sound, the isentropic 
flow requirement, and the equations of state for the 
enthalpy and density we get: 

PoAlV* = PnAuVn. (4) 
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= S m = + x n ,(s^ - s J) (8) 

K = Kj + ** <Kut ~ KJ> (9) 


_1_ = _1_ + x M lj 

^dc PaeJ V Pae,£ ^oe>fj 


(16) 


P nt 


1 

P nt/ 


M 



jJ 

PmJ> 


(10) 



(17) 


The system of equations from Equation (4) through 
Equation (10) can be solved iteratively by guessing the 
pressure at the nozzle throat P* and using Equation (8) 
to compute the quality at the throat x^. The quality can 
then be used to compute the specific enthalpy and 
density at the throat using Equations (9) and (10), 
respectively. Equations (4) and (5) can now be solved 
in the velocities V* and V w . The speed of sound at the 
throat can be computed from Equation (7) and should 
be compared to the velocity term V^. If the two 
quantities are found different in value, another value of 
pressure is tried until convergence. 

At this point, the mass flow rate of the primary 
fluid, hip, may be computed using 


Again, the pressure at the nozzle exit, P^, is assumed 
and the above system is solved by iteration until 
convergence. 

Secondary Flow : 

Since the secondary flow is all liquid, 
Bernoulli’s equation can be applied between the inlet 
and exit of the secondary flow tube: 

P = ♦ - Pib V* (18) 

where it has been assumed that the liquid density is 
nearly constant throughout the secondary fluid path; 


= p,A«v nt <u> 

The procedure described above can now be repeated 
between the throat and exit sections of the primary 
nozzle to produce the quality, density, specific 
enthalpy, Mach number, and velocity at the exit 
section. This can be achieved using the following: 

s = s = s , + x (s -s,) (12) 

nt nt n ej bc^dc^ n ej' 


The continuity equation in the secondary flow 
path yields 

m, = 4»m p = Px A x V k (20) 

where the entrainment ratio, <t\ is known a priori. 
Equations (18) through (20) can thus be solved to 
produce the pressure and velocity at the secondary flow 
path exit, P* and V*, respectively. 


“p = PnAe V « 


K + 



(13) 

(14) 


Mixing Chamber : 

The underlying assumption in the mixing 
chamber is: P^ = P^ = the pressure along the slanted 
sides. The continuity momentum, and energy equations 
in the mixing chamber are thus: 


K = \cj * *n.(V* ' Vf) 


(15) 
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riip (1 + 4>) 



P «= \e + P * A * - *mcKc = 
dip (!♦*) - riipV^ - 


K- - m 

K* + \ y2 j (i + *) 


than unity, a shock wave is assumed to occur and the 
following system of equations is assumed to apply: 

p V = p V (29) 

“ me me r sa u 

p +p vi = P t p (30) 

me • me me n r # as 





(31) 


Equations (21) through (23) along with the equations of 
state 


S n* = 


+ ^me (^me^ ^mej) 

(24) 

K* = 


+ (*w - KJ 

(25) 

1 

1 

.«_( 1 - 1 ) 

(26) 

Pn* 

P mej 

l Pn^ Pn*J 



P - P (PJ 0 ^ 

Again, the pressure is assumed and an iterative solution 
is carried out until convergence is reached. 

If, on the other hand, the Mach number at the 
mixing chamber exit is less than unity, a subsonic 
solution is assumed and an isentropic flow analysis 
similar to the one carried out previously is performed 
in the diffuser. This yields the conditions at the 
diffuser exit, P*, V*, x*. h*, />*, M*, and a*. 


and the equations defining the Mach number and the 
speed of sound 



constitute a system capable of fully determining the 
state of the mixture at the exit of the mixing chamber. 
This includes p„, V TCf A TC , h^. s re , and a^. 

Diffuser : 

The conditions at the mixing chamber exit 
were assumed identical to those at the diffuser inlet for 
purposes of the analysis presented in this paper. Before 
this analysis is performed, a test for the existence of 
shock waves and for whether the flow is subsonic or 
supersonic is carried out. The Mach number computed 
above determines the latter issue and in case it is larger 


RESULTS AND DISCUSSION 

The mathematical model presented in the 
previous section was solved numerically in order to 
determine the distribution of thermodynamic states. In 
addition, the area of the mixing chamber exit was 
calculated as part of the solution instead of being 
specified a priori. This was necessitated by the 
constant-pressure mixing assumption applied to the 
mixing chamber. The analysis of a variable-area, 
fixed-geometry ejector of this type would necessarily 
require a streamwise pressure gradient during the 
mixing process except at a single operating point. In 
other words, the results associated with the constant- 
pressure mixing assumption are a series of design points 
as opposed to an analysis of an ejector of fixed 
geometry. Previous analyses of this type (Fairuzov and 
Bredikhin, 1995) have attempted to analyze fixed- 
geometry ejectors using a quasi constant- pressure 
model. This approach seems to create potential 
contradictions in simultaneously satisfying the constant- 
pressure constraint and the conservation laws in the 
mixing chamber. 
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Figure 2 displays the mixing chamber exit-to- 
diffuser exit area ratio as a function of the entrainment 
ratio for nozzle inlet temperatures of -5, 0, and 5°C. 
As is evident, the area ratio is a rather complicated 
function of the entrainment ratio for specific primary 
and secondary fluid inlet states. This is primarily the 
case because of the competing effects of the 
monotonically increasing velocity and the relatively 
complicated variation in density. For a specified mass 
flow rate, the area varies inversely with the product of 
density and velocity, hence leading to the observed 
trends in Figure 2. The density itself varies due to the 
decreasing pressure in competition with the decreasing 
mixture quality, shown in Figure 3. The latter falls 
rather rapidly with the entrainment ratio in the low 
entrainment ratio regime, while it exhibits a minimum 
at higher values. This behavior is explained by the 
quadratic dependence of the mixing chamber pressure 
on the entrainment ratio. In the low entrainment ratio 
regime, the pressure varies little so that the observed 
reduction in quality is due to the effect of the low 
secondary fluid temperature. On the other hand, for 
higher entrainment ratios, the thermal effect becomes 
dominated by the rapidly decreasing pressure. The 
effect of the primary inlet temperature is thermal, 
similar to the effect of the secondary fluid in the low 
entrainment ratio regime. Thus, lower primary fluid 
temperatures are observed to uniformly produce lower 
quality. 

Figure 4 shows the variation of the mixing 
chamber exit Mach number with respect to the 
entrainment ratio and primary fluid inlet temperature. 
Increasing the entrainment ratio increases the Mach 
number because of the lower mixing chamber pressure 
and the associated higher mixture velocity. While the 
present model is also capable of handling supersonic 
flow cases, in the present paper we present only 
subsonic and sonic cases. The effect of increasing the 
primary fluid temperature, which also implies an 
increase in the primary fluid pressure, is to create a 
more energetic primary jet entering the mixing 
chamber, hence increasing the mixture velocity. An 
interesting observation is seen in the low entrainment 
ratio regime where a minimum mixture velocity 
increases monotonically. In this region, the speed of 
sound increases more rapidly than does the velocity. 
As the primary fluid temperature approaches that of the 
secondary, the primary fluid may not be energetic 
enough to cause the mixture to reach sonic conditions. 
This trend is observed for the -5°C primary fluid 


temperature case. 

The diffuser decelerates the flow from the 
mixing chamber throat to the ejector exit, which is of 
constant area. Shown in. Figure 5 is the variation of the 
diffusers exit quality as a function of the entrainment 
ratio and primary fluid inlet temperature. There are 
competing physical effects of pressure and temperature 
analogous to those described in Figure 3. In the low 
entrainment ratio regime, the thermal effect dominates 
the mild pressure effect, while the reverse is true in the 
high entrainment ratio regime. For all cases examined, 
the ejector exit quality remained above zero in contrast 
to some of the experimental observations reported by 
Fairuzov and Bredikhin (1995) where liquid was 
observed at the ejector exit under certain conditions. 
The constraints of the constant-pressure mixing, 
coupled with the choice of baseline fluid inlet states, 
selected a subset of possible experimental cases, none 
of which corresponded to one of those of Fairuzov and 
Bredikhin (1995). 

Figure 6 shows the variation of the diffuser 
exit Mach number as a function of the entrainment ratio 
and the primary fluid inlet temperature. As entrainment 
ratio is increased, one observes a knee in the curve 
above which the exit Mach number rapidly becomes 
unmanageably high, especially with a low primary fluid 
temperature. In that case, very little diffusion is 
occuring because the throat size necessary to 
accommodate constant-pressure mixing approaches that 
of the diffuser exit. At higher primary fluid 
temperatures, a practical amount of diffusion is possible 
even with choked conditions at the mixing chamber 
exit. 

The variation of the diffuser static and 
stagnation exit pressures as functions of the entrainment 
ratio and the primary fluid inlet temperature can be 
seen in Figures 7 and 8. Similar trends are observed in 
which the pressure ratio decreases with increasing 
entrainment ratio, qualitatively similar to performance 
curves for fixed-geometry ejectors. As expected, for 
the low primary fluid temperature case, the 
compression ratio is greater than one only for very low 
entrainment ratios. On the other hand, the high 
primary fluid temperature case resulted in practical 
compression ratios over the entire range. 
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CONCLUSIONS 

This paper presented a design-oriented model 
for a jet pump employing a two-phase primary fluid 
inducing a secondary liquid having the same chemical 
composition. The model is capable of accounting for 
phase transformations due to compression, expansion, 
and mixing employing an isentropic homogenous 
expansion/compression method of analysis. For given 
system requirements for flow rate and inlet fluid states, 
the model predicts the required mixing chamber throat 
necessary to accomplish constant-pressure mixing. For 
practical cases, the required throat size increases with 
increasing entrainment ratio. Practical compression 
ratios can only be achieved with sufficiently energetic 
primary fluids. Similar to fixed-geometry ejectors, the 
variation of compression ratio with entrainment ratio 
follows a monotonically decreasing trend. The more 
energetic primary fluid appears to be more suitable for 
inducing high compression ratios than high entrainment 
ratios. 
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ni = nozzle inlet nt = nozzle throat ne = nozzle exit 

me = mixing chamber exit de = diffuser exit 


Figure 1. Schematic diagram of ejector with two-phase primary fluid and subcooled liquid 
secondary fluid. 



Figure 2. Mixing chamber throat to diffuser exit area ratio versus entrainment ratio for 
three values of the primary fluid inlet temperature. 
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Figure 5. Diffuser exit quality versus entrainment ratio for three values of the primary fluid 
inlet temperature. 



Figure 6. Diffuser exit Mach number versus entrainment ratio for three values of the 
primary fluid inlet temperature. 
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Figure 7. Ejector compression ratio based on exit static pressure versus entrainment ratio 
for three values of the primary fluid inlet temperature. 





Figure 8. Ejector compression ratio based on exit stagnation pressure versus entrainment 
ratio for three values of the primary fluid inlet temperature. 




